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Q^l Abstract 

, The Coma cluster was observed with an X-ray Imaging Spectrometer (XIS) onboard Suzaku in six 

pointings, including the central X-ray peak region, the 14' west offset region, 30' and 34' north-west 
CN ; offset regions, and 44' and 60' south-west offset regions. Owing to its lower background level, Suzaku has 

better sensitivity to Fe Ka lines than other satellites. Using precise Fe line measurements, we studied 
I the temperature structure, possible bulk motions, and iron abundance distributions in the intracluster 

medium (ICM). The observed spectra were well- represented by a single-temperature model, and two- or 
three- temperature model did not improve substantially. The temperature, derived from Ka line ratios 
. of H-like and He-like Fe, agree with those derived from the single-temperature model. Because the line ratio 

is a steep function of temperature, the consistency supports the accuracy of temperature measurements 
Q , conducted with Suzaku. Within the 34' region, the redshift derived from the central energy of the He-like 

Vh ' Fe line is consistent with that from optical observations, within a calibration error of 18 eV or 818 km s~^ 

C/5 I in the line of sight. This value is smaller than the sound velocity of ICM, which is 1500 km s~^. The 

^ ■ central energy of Fe lines at the 44' offset region around the NGC 4839 subcluster is also consistent with 

those within the 34' region. These results on the temperature and velocity structure suggest that the core 
. of the cluster is in a relaxed state, and non-thermal electrons relevant to the radio halo are accelerated by 

^ ' intracluster turbulence rather than large-scale shocks. Fe abundance is almost constant at 0.4 solar within 

, the 34' region, and decreases with radius. This value is slightly lower than those of other clusters, which 

' means the gas have been mixed well during a past merger associated with the growth of the cluster. 

Key words: galaxies: clusters: individual (Abell 1656, the Coma cluster) — X-rays: galaxies: clusters 
— intergalactie medium 
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1. Introduction effect. Tamura et al. (2011) discovered a significant bulk 

motion in Abell 2256, which is a well known merging clus- 

Clusters of galaxies are thought to grow into larger sys- ter. Schuecker et al. (2004) performed a Fourier analysis 

tems through complex interactions between smaller sys- of XMM- Newton data of the Coma cluster, which revealed 

tems. Signatures of merging events include temperature the presence of a scale-invariant pressure fluctuation rang- 

and density inhomogeneities and bulk motions in the in- ing between 29 and 64 kpc, and '^10% of the total pressure 

tracluster medium (ICM). Numerical simulations predict in turbulent form. However, in relaxed clusters, significant 

simplified bulk motions with a substantial fraction of virial bulk motions and turbulence have not been detected. On 

velocity (~1000 km s^^), lasting several Gyr after each the basis of ASCA and Chandra data, (Dupke & Bregman 

merging event (e.g., Roettiger et al. 1996; Norman & 2001; Dupke & Bregman 2006) claimed a large velocity 

Bryan 1999). Measurements of the ICM temperature and gradient of ^2400 km s~^ in the Centaurus cluster, which 

velocity structure, therefore, provide crucial information is a relaxed cluster with a cool core. However, Ota et al. 

for understanding the evolution of clusters. (2007) found a negative result of |A?;| < 1400 km s~^ on 

Velocity measurements of bulk motion and turbulence the basis of Suzaku observations of the Centaurus cluster, 

are essential for mass estimation and cosmological stud- Sato et al. (2008) also searched for possible bulk motions 

ies. If ICM has a significant bulk velocity compared to in the AWM 7 cluster with Suzaku, and observed no sig- 

its thermal velocity, the associated non-thermal pressure nificant bulk motions. The upper limit of gas velocity was 

would threaten the assumption of hydrostatic equilibrium |Au| < 2000 km s~^. Fujita et al. (2008) do not detect 

in deriving the total gravitational mass in a cluster. For the variation of the redshift of the ICM in the Ophiuchus 

example, factors of 2-3 discrepancies between the X-ray cluster, and the upper limit of the velocity difference is 

and lensing mass estimates in some objects (e.g., Ota et 3000 km s~^. Sugawara et al. (2009) studied the ICM 

al. 2004; Hattori et al. 1999) could be partly due to this flow in Abell 2319 with Suzaku, and detected no veloc- 
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ity differences within the observed region. Using spectra 
with a Reflection Grating Spectrometer (RGS) onboard 
XMM-Newton, Sanders et al. (2011) provided upper lim- 
its on turbulence velocities in cluster core regions: at least 
15 sources had less than 20% of the thermal energy density 
in turbulence. 

The Coma cluster (z = 0.0231), also known as 
Abell 1656, is one of the most studied clusters (Biviano 
1998) and has been observed at all wavelengths from radio 
to hard X-ray bands. However, the physical state, partic- 
ularly the dynamical state, of the Coma cluster has not 
been completely understood. Measurements of the veloc- 
ity dispersion and distribution of galaxies belonging to the 
Coma cluster constrain the dynamical history of the Coma 
cluster. Fitchett & Webster (1987) found that two central 
galaxies, NGC 4889 and NGC 4874 have a significant dif- 
ference in velocity, which is evidence of a recent merger. 
Collcss & Dunn (1996) argued that NGC 4874 was the 
original dominant galaxy of the main cluster, and that 
NGC 4889 belonged to a subgroup, that recently merged 
with the main cluster. A substructure associated with 
the galaxy, NGC 4839, located 40' south-west of the clus- 
ter core, was found by Mellier et al. (1988) and Merritt 
& Trimbley (1994). The existence of this subcluster was 
confirmed by Collcss & Dunn (1996). 

The thermal X-ray emission of the Coma cluster has 
been observed with several X-ray satellites. Briel et al. 
(1992) constructed an X-ray map of the Coma cluster 
using ROSAT All- Sky- Survey data and determined the 
X-ray surface brightness profile of the cluster out to a ra- 
dius of roughly 100'. The temperature maps of clusters 
provide knowledge about the history of past subcluster 
mergers. On the basis of ASCA observations, Honda et 
al. (1996) and Watanabe et al. (1999) found that ICM is 
not isothermal. Arnaud et al. (2001) used XMM-Newton 
to study the temperature structure in the central region of 
the Coma cluster. The projected temperature distribution 
around NGC 4889 and NGC 4874 is remarkably homoge- 
neous, which suggests that the core is mostly in a relaxed 
state. Except at the center, the temperature decreases 
slightly with radius. A cool filament of an X-ray emission 
in the direction of the galaxy NGC 4911, which is located 
in the south-east of the cluster center, was detected with 
Chandra (Vikhlinin et al. 1997) and XMM (Arnaud et 
al. 2001). ROSAT observed a substructure around the 
NGC 4839 subgroup (White et al. 1993), and a hot region 
in the direction of NGC 4839 was observed by Arnaud 
et al. (2001). By XMM-Newton observations, Neumann 
et al. (2001) found compelling evidence for the subgroup 
around NGC 4839 was on its first infall into the Coma 
cluster, but it had not passed its core yet. This interpre- 
tation is different from that of Burns et al. (1994), who 
suggested that the NGC 4839 group was moving out of 
the cluster after having already passed through the Coma 
cluster. 

Non-thermal electrons, observed via diffused radio syn- 
chrotron emission, have been detected in more than 50 
clusters, and all of them are undergoing mergers (e.g., 
Buote 2001; Schuecker et al. 2001). The Coma cluster 



also has a cluster- wide synchrotron radio halo, emitted by 
relativistic electrons due to merging (Feretti & Giovannini 
1998). Recent non-thermal detections have been claimed 
by Rcphaefi & Grubcr (2002) with RXTE and by Fusco- 
Femiano et al. (2004); Fusco-Femiano et al. (1999) with 
BeppoSAX, although the latter detection is controversial 
(Rosetti & Molendi 2004; Fusco-Femiano et al. 2007). 
Long observations (~1 Ms) by INTEGRAL have imaged 
extended diffuse hard X-ray emission from the Coma clus- 
ter, although it was found to be completely consistent 
with thermal emission (Renaud et al. 2006; Eckert et al. 
2007; Lutovinov et al. 2008). Using the temperature map 
of ICM obtained with the FN detector onboard XMM- 
Newton, Wik et al. (2009) analyzed data obtained from 
the Suzaku Hard X-ray Detector (HXD) (Takahashi et 
al. 2007) and derived the strongest upper limit of non- 
thermal emission. 

This paper reports results from six Suzaku observations 
of the Coma cluster out to 60', ~1.5 Mpc, conducted using 
the X-ray Imaging Spectrometer (XIS) (Koyama et al. 
2007) onboard Suzaku (Mitsuda et al. 2007). Owing to 
its lower background level, Suzaku has better sensitivity 
to the iron Ka lines than previous satellites. In addition, 
the accurate calibration of XIS allows us to measure the 
velocity of ICM precisely. 

We used the Hubble constant, Hq — 70 km s~^ Mpc~^. 
The distance to the Coma cluster is Dl = 101 Mpc, and 1' 
corresponds to 28.9 kpc. We used the abundance ratio by 
Lodders (2003), in which the solar Fe abundance relative 
to H is 2.95xl0~^. Errors are quoted at a 90% confidence 
level for a single parameter. 

The remainder of this paper is organized as follows: we 
present the observations in section 2 and describe the data 
analysis in section 3. We describe the results of spectral 
fittings and temperature structure in subsection 4.1, bulk 
motions in subsection 4.5, and the Fe abundance in sub- 
section 4.6. In section 5, we discuss the results, and in 
section 6, we summarize the findings. 

2. Observation and Data Reduction 

Suzaku carried out six observations in the Coma cluster. 
The details of the observations are summarized in table 1. 
Two central observations were carried out in May 2006, 
during the Suzaku Phase-I period. The first observation, 
"Coma Radio halo" , has a pointing that is 5' offset from 
the X-ray peak of the Coma cluster, which has coordinates 
(RA, Dec) = (194.9367, 27.9472) in degrees. The second 
observation, "Coma cluster offset", was located 14' west 
offset from the first observation. The other observations, 
in the 30' and 34' north west offset regions and 44' and 
60' south west offset regions, were carried out in 2007, 
during the Suzaku Phase-II period. To eliminate back- 
ground emissions, archival Suzaku data with a 5 degree 
offset from the Coma cluster, observed in June 2007 was 
also used, and is shown in table 1. 

XIS was operated in its normal mode during observa- 
tions. The XIS instrument consists of four sets of X-ray 
CCD (XISO, 1, 2, and 3). XISl is a back-illuminated (BI) 
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Table 1. Suzaku observation log of the Coma cluster 



Field name 


Target name 


Sequence 


Observation 


RA 


Dec 


Exposure 


Distance from the 






number 


date 


(dcg) 


(dcg) 


(ks) 


X-ray peak* 


center 


Coma Radio Halo 


801097010 


2006-05-31 


194.9267 


27.9061 


150 


4.'7 


14' offset 


Coma cluster offset 


801044010 


2006-05-30 


194.6939 


27.9466 


79 


13'.9 


30' offset 


Coma 11 


802082010 


2007-06-19 


194.6305 


28.3939 


53 


30.'1 


34' offset 


Coma BKG2 


802084010 


2007-06-21 


194.3428 


28.1403 


30 


33.'7 


44' offset 


Coma 45 


802047010 


2007-12-02 


194.2558 


27.5714 


30 


44f4 


60' offset 


Coma 60 


802048010 


2007-12-04 


194.0251 


27.4252 


35 


59f5 


5° offset 


Coma BKG 


802083010 


2007-06-21 


198.7472 


31.6480 


30 


4.°9 



* X-ray peak has coordinates (RA, Dec) = (194.9367, 27.9472) in degrees. 





Fig. 1. XISO image (0.5—4.0 keV) of the Coma cluster. Left: Differences in exposure time and the vignetting effect are uncorrected. 
Regions of spectral accumulation are shown as blue, red, and green squares. The light blue circles show the locations of NGC 4889, 
NGC 4874, and NGC 4839. Right: Definition of the cell numbers in the center and 14' offset regions. The cell numbers 12 and 15 
of the center region (green numbers) were excluded from our analysis because of strong emissions from calibration sources. 



sensor, while XISO, 2, and 3 are front-illuminated (FI) 
sensors. XIS2 was not used during the four offset obser- 
vations beyond 30' from the X-ray peak. Figure 1 shows 
a 0.5-4.0 keV image of XISO. 

We performed data reduction using HEAsoft version 
6.6.3. XIS event lists created by the rev 2.0 pipeline pro- 
cessing were filtered using the following additional criteria: 
the geomagnetic cut-off rigidity (C0R2) > 6 GV, and ele- 
vation angle from the earth limb > 10°. The data formats 
of 5x5 and 3x3 editing modes were added. The exposure 
times after the data selection are shown in table 1. 

To study the temperature structure and search for pos- 
sible bulk motions of ICM in the Coma cluster on the scale 
of a few arcminutes, we divided the 18' x 18' square XIS 
field of view (FOV) of the center and 14' offset regions 
into 4x4 and 2x2 cells, respectively, as shown in figure 
1. Each region in the center and 14' offset regions sub- 
tends 4f5x4.'5 (117x117 kpc), and 9' x 9' (234x234 kpc), 
respectively. Regions around calibration sources were ex- 
cluded. Luminous point sources in which the flux limit 



corresponds to ^ 10~^ counts s~^ in the 2.0-10.0 keV en- 
ergy range were also excluded as circular regions with a 
radius of 1'. At most 10 point sources were excluded for 
each of the observed regions. We extracted spectra from 
the whole XIS FOV for regions more than > 30' from 
the X-ray peak. The non X-ray background (NXB) was 
subtracted from each spectrum using a database of night 
Earth observations with the same detector area and COR 
distribution (Tawa ct al. 2008). 

We included the degradation of energy resolution due to 
radiation damage in the redistribution matrix file (RMF) 
generated by the xisrmf gen Ftools task. We also created 
an ancillary response file (ARE) using the xissimarf gen 
Ftools task (Ishisaki et al. 2007). A decrease in the 
low-energy transmission of the XIS optical blocking filter 
(OBF) was also included in ARE. To generate ARE files 
for the center and 14' offset regions, we used a /3-model 
profile for the simulated surface brightness profile in Briel 
ct al. (1992). For 30', 34', 44', 60', and 5 degree offset re- 
gions, we generated AREs assuming a uniformly extended 
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emission from an encircled region with a 20' radius, be- 
cause the gradient of surface brightness within each FOV 
is smaU. We used the XSPEC_vll.3.2ag package for spec- 
tral analysis. 

3. Data analysis 

We have investigated the temperature structure, bulk 
motions, and iron abundance distribution of the Coma 
cluster. We fitted the spectra with a thermal plasma 
model (APEC: Smith et al. 2001) to obtain the tempera- 
ture, rcdshift and Fc abundance of ICM. A detailed spec- 
tral analysis is presented in subsection 3.2. We used three 
Gaussian models to obtain the temperature using the nor- 
malization ratio of Ka lines of H-likc and He-like Fc, and 
rcdshift using the shift of line centroid energy, which is 
presented in subsection 3.3. 

3.1. Estimation of background spectra 





XISO 
XIS1 

XIS3 

1 1 1 — 1 — 







_i I i_i : 1 1 1 L. 

0.5 1 2 5 

Energy [keV] 



Fig. 2. NXB-subtracted spectra of the 5 degree offset region 
in the 0.3-8.0 (red: XISl) and 0.5-8.0 keV energy ranges 
(black and green for XISO and 3, respectively), ignoring the 
1.40-1.55 keV energy range, where a relatively large uncer- 
tainty exists in the instrument background (Tawa et al. 2008). 
The bottom panel shows the residuals of the fit. 

First, we fitted the spectra in the 5° offset region to 
determine the local X-ray background. As shown in 
Yoshino et al. (2009), the background emission of Suzaku 
XIS can be fitted with a three component model: two 
thermal plasma models (APEC: Smith et al. 2001) for 
the local hot bubble (LHB) and the Milky Way halo 
(MWH) , and a power- law model for the extragalactic cos- 
mic X-ray background (CXB). MWH and CXB compo- 
nents were convolved with an absorption in the Galaxy. 
We, therefore, fitted the spectra simultaneously in the 0.5- 
8.0 keV energy range for XIS FIs and the 0.35-8.0 keV 
energy range for XIS BI using the following model for- 
mula: apeci,}iB +W0'bs * (apecMWH+power-lawcxB)- We 
assumed a zero redshift and a solar abundance for LHB 
and MWH components. The temperature of LHB was 
fixed at 0.07 keV (Takei et al. 2008), while the temper- 
ature of MWH was variable. The column density of the 
Galactic neutral hydrogen was fixed at 1.0x10^° cm~^ 



(Kalberla et al. 2005). To avoid systematic uncertainties 
in the background, we ignored the 1.40-1.55 keV energy 
range, where a relatively large uncertainty exists because 
of an instrument Al line (Tawa et al. 2008), and energies 
above 8 keV. We also excluded the narrow energy band 
between 1.82 and 1.84 keV in the fits because of incom- 
plete responses around the Si edge. Results of the spectral 
fit are shown in figure 2 and the resultant parameters of 
the fit are shown in table 2. The derived photon index and 
normalization of the CXB component agree with Kushino 
et al. (2002) and Takei et al. (2008). The derived temper- 
ature and normalizations of the two thermal models are 
consistent with those obtained from XMM and Suzaku 
observations (Lumb et al. 2002; Yoshino et al. 2009). 

We generated the simulated spectra of the X-ray back- 
ground including LHB, MWH, and CXB emissions derived 
from the resultant parameters of the fits in the 5° offset 
region. Wc used the simulated spectra as the background 
to fit the spectra for ICM of the Coma cluster. Figure 
3 shows the comparison between the XISO spectra of the 
60' offset region of the Coma cluster and backgrounds. 
Above 7 keV, NXB dominates the observed spectra of the 
60' offset region. We also show the background-subtracted 
spectrum of the 60' offset region of the Coma cluster in 
figure 3 (blue). The X-ray background spectra are well 
reproduced. 

3.2. Spectral fits with single-, two- and three temperature 
models 

We fitted the XIS spectra in each region of the Coma 
cluster using a single temperature model (APEC) with 
the Galactic absorption, Nn. Each spectral bin contained 

50 or more counts. We fitted the spectra in two energy 
bands: 1.0-8.0 and 5.0-8.0 keV. To avoid systematic un- 
certainties in the background, we ignored energies higher 
than 8 keV. We also excluded the narrow energy band 
between 1.82 and 1.84 kcV in the fits because of incom- 
plete responses around the Si edge. Temperature, redshift 
and normalization of the single temperature model were 
free parameters, and A'h was fixed to the Galactic value, 
1.0 X 10^° cm~^, in the direction of the Coma cluster. 
Abundances of He, C, N, and Al were fixed to be a solar. 
We divided the other metals into three groups: O, Nc and 
Mg; Si, S, Ar and Ca; and Fe and Ni, and allowed them 
to vary. The resultant parameters are summarized in ta- 
ble 3, and a sample of the spectra is shown in figure 4. 
Since our data have very high statistics, especially in the 
central region, there remain residual structures around the 

51 edge structure. These residuals are probably caused 
by systematic uncertainties in the response matrix be- 
cause the BI and FI detectors gave different residuals. 
Therefore, to avoid systematic uncertainties in all regions, 
we excluded the energy range of 1.7-2.3 keV and fitted the 
spectra again. The derived temperatures hardly changed, 
although the reduced decreased by 20%. The resultant 

are summarized in table 3. 

To constrain the temperature structure, wc fitted the 
spectra of the center and offset regions with a two- 
temperature (APEC) model in the energy range of 1.0-8.0 
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Cosmic X-ray Background Local Hot Bubble Milkyway Halo 

r normalization* Temperature normalization^^ Temperature normalization^^ Reduced-x'^ 

(kcV) (keV) xVd.o.f^ xVd.o.f§ 

1.38±0.07 7 6±0.3 0.07 (fixed) 6.79±2.57 0.16±0.03 6.55±2.57 188/141 223/154 

* Units: photons cm~"^ s~^ keV^^ sr^^ at 1 keV. 

t Normalization of the APEC, component divided by the solid angle, fi", assumed in the uniform-sky ARF calculation (20' radius), 
Norm = ^ nenjjdV / + z)^ D'^) /Q,^ X 10~^* cm~^ (4007r)~^ arcmin"^, where Da is the angular distance to the source, 
-f Resultant x^/d.o.f fitted excluding 1.40-1.55 keV. 
§ Resultant x^/d.o.f fitted including 1.40-1.55 keV. 



60' offset 



> 

0) 



o 
O 




Energy [keV] 

Fig. 3. XISO spectra of the 60' offset region of the Coma cluster (black), NXB (green), and the simulated background spectrum 
including LHB, MWH and CXB (red). Blue crosses show the background-subtracted spectrum of the 60' offset region of the Coma 
cluster. 



keV, ignoring the energy range of 1.7-2.3 keV. The are 
improved slightly by less than a few percent. The derived 
temperatures of one component are similar to those de- 
rived from the single-temperature model fits, and those 
of the other component are either above 10 keV or below 
2 keV. To obtain more detailed constraints for the lower 
and higher temperature components, we fitted the spec- 
tra with a three-temperature (APEC) model in the energy 
range of 1.0-8.0 keV, excluding the 1.7-2.3 keV energy 
range to avoid uncertainties in the response matrix. In 
this model, the lowest, middle and highest temperatures 
are restricted to below 6 keV, within 6.0-10.0 keV, and 
above 10.0 keV, respectively. The resultant values are 
shown in table 3. 

To derive the Fe abundance and velocity of the bulk- 
motion from Fe-K lines, we fitted the spectra of each re- 
gion with the single-temperature APEC model in an en- 
ergy range of 5.0-8.0 keV. The results are shown in table 
3. 

To precisely determine the bulk velocity of ICM, the 
accuracy of photon energy measurements is crucial. Some 
sensors have large uncertainties in the determination of 
the line centroid energy. Thus, we averaged the spectrum 
over FOV in (i) all detectors, and (ii) all detectors except 
for different Mn Ka line centroid energies to check the 
spectral shape and measured redshift in the center and 



NGC 4839 regions (44' offset region). 

3.3. Central energy and normalization ratio oj Fe lines 

To derive the ICM temperature from the normalization ra- 
tio of Kof lines of H-like and He-like Fe, and the possible bulk 
motions of ICM from the central energy of the lines, we fit- 
ted the spectra in the energy range of 5.0-8.0 keV with a 
sum of bremsstrahlung and three Gaussian components, which 
correspond to Kq lines of He-like and H-like Fe, and a mix- 
ture of the Ka line of He-like Ni and the K/3 line of He-like 
Fe. The temperature and normalization of the bremsstrahlung 
component, and the central energy and normalizations of the 
Gaussian components were free parameters, while the line 
width of Gaussian models was fixed at 0. The best-fit spectra 
are shown in figure 5. 

4. Results 

Temperature structure of ICM 

Figure 6 shows a radial profile and map of the derived tem- 
peratures from the spectral fits with the single temperature 
(APEC) model in the energy range of 1.0-8.0 keV. The ra- 
dial temperature profile of the Coma cluster is relatively flat 
within ~500 kpc, which corresponds to the 14' offset region, 
and decreases with the radius to the outer region. While the 
regions observed with Suzaku are limited, the resultant tem- 
perature profile is similar to that derived from XMM-Newton 
observations, which cover nearly the entire Coma cluster out 
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center (cell numbers) 14' offset (cell number 17) 




Energy [keV] Energy [keV] 

Fig. 4. Spectra for one cell in the center and 14' offset regions, and all otfier offset regions, fitted by the single temperature (APEC) 
model. Black, red, blue, and green crosses and lines correspond to the XISO, 1, 2, and 3, respectively. The middle and bottom panels 
show the residuals of the fit. 
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Fig. 5. Representative iron Ka line spectra (for the same regions as in figure 4), fitted with the bremsstrahlung and three Gaussian 
model in the energy range of 5.0—8.0 keV. Three Gaussian models (black lines) represent the strong Ka line of He-like Fc, a weaker 
line of H-like Fc, and a mixture of the Ko line of He-like Ni and the K/3 line of He-like Fe. 
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Table 3. Resultant parameters of the fits with the single-temperature (APEC) model in the energy ranges of 1.0-8.0 keV and 
5.0-8.0 keV, the two- or three-temperature model, and three Gaussians models. In the case of three Gaussian models, we converted 
the normalization ratio of Ka lines of H-like and He-like Fe to temperature with the APEC model. 







APEC fit in 1.0-: 


B.O keV 




APEC fit in 5.0-5 


5.0 kcV Gaussians 


Region 


kT 
[kcV] 


xVd.o.f.* xVd.o.f.T 
(IT) (IT) 


xVd.o.f.^ 
(2T) 


" xVd.o.f.i* 

(3T) 


Fe Redshift 
[solar] [10-2] 


X^/d.o.f. kT from line 
ratio [keV] 



center 



8.47l:°o8 938/829 

1 7.8lto'?o 1373/1108 

2 7.19t[!J^ 1325/1108 

3 7.12;°'}^ 901/829 

4 8.73t[j'ie 1363/1108 

5 8.20j:[j;]^ 1574/1108 

6 7.8i;^'°^ 1298/1108 

7 7.45t[!;}2 1314/1108 

8 8.93;^!^ 1290/1108 

9 8.66+o'?2 1533/1108 

10 8.26tg'!^^ 1395/1108 

11 8.15;^'ii^ 1379/1108 

13 9.08l:^'23 1 361/1108 

14 8.95lo'i^ 1306/1108 
14'ofrset 





16 
17 
18 
19 


9.141°:!^ 

9 22+°-i^ 

9.oo;°i^ 

8.43tg:l? 


1355/1108 
1201/1108 
1303/1108 
1280/1108 


1047/950 
987/950 
1087/950 
1067/950 


1035/948 
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to ~1200 kpc (Matsushita 2011). 

Excluding the south-east quadrant (cell numbers: 2,3,6, and 
7), and the 14' offset region, temperatures of cells in the cen- 
ter region have a flat distribution within a small range between 
~8 keV and ~9 keV. Within 300 kpc from the X-ray peak in 
the center region, the weighted mean temperature with statis- 
tical errors of the 16 cells is 8.34±0.19 keV. The south-east 
quadrant (cell numbers: 2,3,6, and 7) in the center region has a 
slightly cooler temperature: the weighted average with statisti- 
cal errors of the four cells is 7.39±0.11 keV. This value is close 
to that of 5.0-7.0 keV in the cool temperature region observed 
by XMM-Newton, which would be a counterpart of the south- 
east quadrant region (Arnaud et al. 2001). On the other hand, 
the north-west quadrant in the center region (cell numbers: 
8,13, and 14) and the nearby area of the 14' offset region (cell 
numbers: 17,18) have slightly higher temperatures of ~9 keV. 
The 34' offset region located to the north-west of the 14' offset 
region, also has a relatively high temperature of ~8-9 keV. 
This direction corresponds to the hot spot, kT = 12. 7+3-5 ^'^^ 
in Donnelly et al. (1999) with ASCA, and fcr = 8.4±0.4 keV in 
Arnaud et al. (2001) with XMM-Newton. Recent INTEGRAL 
observations have shown a surface brightness excess relative to 
XMM observations, which is well-represented by the extended 
hot thermal emission with ~ 12 ± 2 keV (Eckert et al. 2007). 



The temperature in the 30' offset region is around ~6-7 keV, 
and is clearly cooler than temperatures in the inner region. 
This value is consistent with the previous results in Honda et 
al. (1996), and Takei et al. (2008). 

The 44' and 60' offset regions include the NGC 4839 sub- 
cluster, and temperatures decreases with radius to 5.31±0.20 
keV and 3.74±0.13 keV, respectively. These values are cooler 
than those in Honda et al. (1996) with ASCA, but are con- 
sistent with those in Wik et al. (2009) and Neumann et al. 
2001 with XMM-Newton. Temperatures of the core and tail 
of NGC 4839 with XMM-Newton, which is located in the 44' 
offset region with Suzaku, are 3.112 4 keV and 4.8I4-Q keV, 
respectively (Neumann et al. 2001). The temperature of the 
main body of the subcluster is 4.4+4-y keV (Neumann et al. 
2001). Suzaku provided more significant temperatures than 
previous measurements in these regions. 

In summary, the temperature distribution observed with 
Suzaku is fairly consistent with the previous results observed 
with ASCA, XMM and INTEGRAL. In addition, we note that 
statistical and systematic errors with Suzaku are much smaller 
than those in previous results. 
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Fig. 6. Radial temperature profile from the X-ray peak (left) and temperature map (right) from our observations of the Coma 
cluster. Temperatures are derived from the spectral fits with the single temperature (APEC) model in the energy range of 1.0-8.0 
keV. The color notations in the left panel are same as those in figure 1. Clear and shaded circles correspond to the center and 14' 
offset regions, respectively, and gray diamonds indicate the results from XMM-Newton (Matsushita 2011). The magenta circle in 
the right panel shows the X-ray peak of the Coma cluster. The temperature map is overlaid on X-ray contours on a linear scale in 
the 0.4-2.4 kcV energy range from the ROSAT- All-Sky-Survey. 




of Single-temperature model 



Fig. 7. Left: Comparison of fitted with single and three temperature (APEC) models. Color notations are the same as those 
in figure 1. Right: Residuals of the fit of XISO with single-temperature (black) and three-temperature (red) models versus energy. 
Each panel shows the region of the spectral fit. 
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4-.2. Fits with three-temperature models 

As shown in figure 4, the spectra in the energy range of 1.0- 
8.0 keV are well-represented by the single-temperature model. 
Except for discrepancies between FI and BI detectors around 
the Si-edge, there are no systematic residuals. Figure 7 shows 
a comparison of the resultant fitted with single and three 
temperature (APEC) models in the energy range of 1.0-8.0 
keV, excluding the 1.7-2.3 keV energy range. A more statisti- 
cally rigorous approach would be to apply the f-test to compare 
single-temperature and three temperature models. The f-test 
indicates that the three-temperature model is significantly bet- 
ter than the single-temperature model in regions 1, 5, 6, 7, 9, 
10, and 11 of the center region and 17 of the 14' offset region, 
and better in regions 2, and 14 of the center region, and 16, and 
19 of the 14' offset region, 30' offset region, 44' offset region, 
and 60' offset region. The largest f statistic value is 11.4 in cell 
number 7 of the center region, and the least f statistic value 
is 0.238 in cell number 8 of the center region. However, the 
single-temperature model fit still represents the spectra fairly 
well. 

4-.3. Temperatures derived from Fe line ratios 

The spectra accumulated over each region may contain 
multi temperature components, although the single temper- 
ature (APEC) model represented those spectra well. Because 
the ratio of He-like and H-like Fe Ka lines strongly depends on 
the plasma temperature, comparisons of temperatures derived 
from the line ratio and spectral fits with the thermal model are 
useful for understanding the temperature structure of ICM. 
Figure 8 shows the ratios of the resultant normalizations of 
the Fe Ka lines, plotted versus the temperature derived from 
the spectral fits with the single temperature (APEC) model 
in the energy range of 1.0-8.0 keV. To convert line ratios to 
plasma temperature, we generated mock spectra assuming the 
single temperature APEC or MEKAL (Mewe et al. 1985; Mewe 
et al. 1986; Kaastra 1992; Liedahl et al. 1995) plasma code 
models, convolved with XIS energy resolution. We then fitted 
the mock spectra with a sum of a bremsstrahlung and three 
Gaussian models. The theoretical line ratios from APEC and 
MEKAL plasma codes are also plotted in figure 8. The de- 
rived line ratios agreed with these theoretical line ratios of the 
APEC plasma code assuming the single temperature model, 
while the MEKAL model gave a line ratio value several % 
lower at a given temperature. Assuming a single temperature 
plasma, we converted the observed line ratios to plasma tem- 
peratures using the theoretical relation for the APEC model. 
Results are shown in table 3 and figure 8. As shown in figure 8, 
line temperatures agreed with those derived from the spectral 
fits very well. 

4-.4- Temperatures at the edge of the radio halo 

Brown & Rudnick (2011) claimed a significant temperature 
difference of ICM at the edge of the radio halo on the basis 
of XMM-Newton observations: The temperature at the inner 
edge of the radio halo is 6.8±0.6 keV, while the temperature 
at the outer edge of the radio halo is 16.3±2.9 keV. To ex- 
amine this temperature difference, we extracted spectra from 
counterpart regions in the 34' offset region, which are repre- 
sented by two red boxes in figure 9, and fitted the spectra with 
the single temperature (APEC) model. As shown in figure 9, 
there are no significant temperature differences for each box in 
the 34' offset region. We then divided FOV of the 34' offset 
region into 2x2 cells as shown in figure 9, and fitted the spec- 



tra with the single temperature (APEC) model. The derived 
temperatures are also consistent with each other. 

4.5. Search for ICM Bulk Motions 

4-5.1. Systematic uncertainty in the energy scale 

To constrain the line-of-sight gas motion using the He-like 
Fe-Ka line, precise calibration of the XIS energy scale is cru- 
cial. Because the six data sets used in the present analysis 
have three different observation epochs (May 2006, June 2009 
and December 2009), it is necessary to check the energy scale 
for each observation. In addition, the energy gain may vary 
from place to place on the same CCD chip owing to the charge 
transfer inefficiency. We thus estimate the XIS calibration un- 
certainty in the following two ways: (i) measurements of the 
calibration source energy and (ii) comparisons of redshift val- 
ues derived from four or three sets of XIS. 

(i) The fiducial absolute energy scale can be examined by 
measuring the centroid energy of the Mn Ka line from built-in 
calibration sources, which illuminate two corners of each XIS 
chip. The spectra accumulated from the calibration source re- 
gions were fitted to the power-law model for continuum and 
two Gaussian functions for the Mn Ka line at 5.894 keV and 
Mn K/3 line at 6.490 keV. Because the Mn K/? line has an 
intrinsically complex line shape and is contaminated by the 
cluster iron emission, only the Mn Ka line is used for the cal- 
ibration. Figures 10 and 11 show the result of spectral fitting 
and the Mn Ka line centroid energy, respectively. 

As shown in the left panel of figure 11, measured line energies 
for the center or 14' offset regions are systematically lower than 
the Mn Ka line energy by 7.3 eV if they are averaged over 
four sensors and two detector segments in each pointing. For 
XIS3, centroid energies are lower than those of the other three 
chips, however, we confirmed that the spectral analysis of the 
redshift measurement with and without XIS3 gave statistically 
consistent results (see section 4.5.3 for more details). We then 
estimate the la systematic uncertainty on the energy scale to 
be (Jsys,Mn = 7 eV. This is justified also from the fact that 
if we calculate the value using the measured Mn Ka line 
centroid Ei and the Icr statistical error Ui in the XlS-i detector, 
and the Mn Ka line energy in literature {E). is defined as 

= Escg Elo(^' - (E))V('^? + <s,Mn). With a,y,,M„ = 

7 eV, becomes less than 10 for 7 degrees of freedom, which 
means the probability of exceeding versus the number of 
degrees of freedom becomes less than 10%: the number of bins 
is 8 for four sensors and two detector segments, and E can 
be regarded as being consistent with (i?) when we consider 
Csys.Mn = 7 eV. In the same manner, the systematic errors are 
obtained as crsys,Mn = 1 eV, 1 eV, 6 eV, 4 eV for 30', 34', 44', 
and 60' offset regions, respectively. Therefore, for the safety's 
sake, we assign 7 eV for the la (68% confidence level) gain 
uncertainty. 

(ii) Though calibration sources provide information on the 
absolute energy scale at the corners of CCD chips, the gain 
may be dependent on position. This intrachip variation can 
be effectively studied by comparing line energies of the same 
sky regions (which correspond to different detector regions) on 
four or three XIS sensors (see also section 3.1 in Ota et al. 
2007). To examine this issue, we divided the center region into 
16 cells, and fitted the APEC model in the energy range of 
5.0-8.0 keV for each XIS sensor. We then estimate the la sys- 
tematic error with a reference to obtained redshift, calculating 

asx' = E^loELo(-f^.(-^--')-(^-.))V(^f'^-"^' + '^s%)for 

each pair of XIS sensors, where E (xis-i) is the measured line 
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Temperature obtained APEC with 1-8 keV 



Fig. 8. Left: Comparison between ratios of tiie normalizations of Fe Ko lines of H-lilie and He-like Fo and the temperatures derived 
from the spectral fits with the single temperature (APEC) model in the energy range of 1.0-8.0 keV. Closed blue and red circles 
with solid error bars correspond to the results of the north and south cells in the center region, respectively. Those with dashed error 
bars show the results of northeast and southwest cells in the 14' offset region, respectively. Green circles indicate the results of the 
four offset regions. Solid and dashed black curves show the theoretical relation between the line ratio and plasma temperature from 
APEC and MEKAL codes, respectively. Right: Comparison of temperatures derived from line ratios of Fe Ka lines and spectral fits 
with the single temperature (APEC) model in the energy range of 1.0-8.0 koV. 




centroid for the XlS-i detector in each region, {E^^) is the 
mean of the line centroid for aU XlS-i detectors in each region, 
and 0"^^^^"'^ is the Icr statistical error in the XlS-i detector. 
With CTsys =11 eV, x'^ becomes less than 17 for 11 degrees of 
freedom, which means the probability of exceeding ver- 
sus the number of degrees of freedom becomes less than 10%: 
the number of bins is 12 for the spectral regions except for the 
calibration source region in the center region. <y\^'^ is ~ 
10 eV, which is larger than the gain uncertainty: thus, we ig- 
nored to avoid underestimating the systematic error. 
Because cr^'^^'^ '' is ignored in these estimates, the obtained 



value is an upper limit for intrachip variation within a 68% 
confidence level. 

In summary, the intrachip variation is an upper limit and is 
larger than the gain uncertainty. Thus, we concluded from (i) 
and (ii) that the Icr systematic error within a 68% confidence 
level is 11 eV. The systematic error is then 18 eV, which cor- 
responds to 818 km s~^ in the line of sight velocity, when we 
quote a 90% confidence level. 
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Fig. 10. Observed XISO spectrum (black crosses) in the energy range around the calibration source for the center (left) and 30' 
offset (right) regions fitted with power-law and two Gaussian models. The red line shows the best fit model and the dashed and 
solid gray lines indicate the bremsstrahlung and Gaussian models, respectively. Mn-Ka, Mn-K/3 lines, and Fe lines from the cluster 
emission are seen in this spectrum. The lower panels show the residual of the fits. 
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Fig. 11. Left: Resultant residuals of the central energy of the Mn-Ko line (5.894 keV) for the center and 14' offset regions. Clear 
and shaded symbols show the segments A and D on XIS, respectively. Right: Same as the left panel for the 30', 34', 44', and 60' 
offset regions. 



4-5.2. Redshijt Measurements from Two Different 
Spectral Models 

As mentioned in subsection 3.2 and 3.3, the line-of-sight ve- 
locity of ICM or redshift can be evaluated using two different 
spectral models: the APEC model or the bremsstrahlung with 
three Gaussian lines. It is then useful to compare these results 
and see if there is any systematic uncertainty in the redshift 
measurement due to spectral modeling. In the APEC model, 
the temperature dependence of the Fe Kq line energy is already 
implemented in the plasma code, which in turn evokes a pa- 
rameter coupling between the temperature and line energy (or 
redshift). On the other hand, in the bremsstrahlung-|-Gaussian 
model, the line centroid should be determined almost indepen- 
dently of the continuum temperature, although its temperature 
dependence has to be corrected later. We did this correction 
by deriving the rest-frame Fe Ka centroid energy as a function 
of temperature on the basis of spectral simulations, assuming 
the APEC model and XIS energy responses. 

In figure 12, the comparison of redshift measured by the 
two different models is shown for all analysis regions defined 
in figure 1. Because there consistency between the two within 
their statistical errors, we conclude that either method ro- 
bustly determines the ICM redshift for the present observa- 
tions. Therefore, we will show results of the fits with the APEC 
model because the uncertainties with APEC are smaller. 



4.5.3. Redshift Measurement from the APEC model 

The obtained redshifts from the spectral fits with the sin- 
gle temperature (APEC) model in the energy range of 5.0-8.0 
keV are shown in figure 12, where the ICM velocity in the line 
of sight is also calculated as vi = c{z — Zd). The resultant red- 
shifts are consistent with the optical redshift in Colless & Dunn 
(1996) from the NED data base^, within statistical and system- 
atic errors, which corresponds to ±0.0027 in the redshift. The 
direct measurements of the central energy of the Ka line of 
He-like Fe also give almost the same redshifts derived from the 
spectral fits with the single temperature (APEC) model. To 
quantitatively derive the redshift, we calculated the mean red- 
shift of the ICM in aU pointings to be (zots) = 0.0227± 0.0005. 
This is slightly lower than the optical redshift, = 0.0231 from 
the NED database, but within 90% error. As shown in table 3, 
in the center and 14' offset regions, the derived highest and low- 
est redshift are Zhighost = 0.0242±0.001 for cell numbers 5 and 8 
of the center region, and ziowest = 0.0209±0.001 for cell number 
2 of the center region. Based on these redshift measurements, 
the derived ICM bulk velocities are Whighcst = 330 ± 300 km 
s~^, and wioweat = —660 ±300 km s~^. Including systematic er- 
rors of ±818 km s~^ with 90% confidence level (Section 4.5.1), 
there is no significant difference. We estimate the 90% upper 

^ http:/ /nedwww. ipac.caltech.edu/ 
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Fig. 12. Left: Comparison of redstiifts derived from spectral fits with tlie single temperature (APEC) model in the 5.0—8.0 keV 
energy range, and the shifted centroid energy of the He-like Fe K« line. Right: Derived redshifts with the single temperature (APEC) 
model in the 5.0-8.0 keV energy range, plotted versus the distance from the X-ray peak of the Coma cluster. Color notations are 
the same as in figure 1. The solid line corresponds to the optical redshift of the Coma cluster, and the dotted and dash-dotted lines 
show the redshift of ±1000 and ±2000 km, respectively. 



limit on the velocity gradient to be Aw = 990 ± 424 ± 1157 km 
s~^. If we further add systematic and statistical errors in the 
quadrature, the upper limit of the gas velocity within the 90% 
confidence level becomes |At;| < 2200 km s~^. For simplicity, 
assuming that the gas is rigidly rotating at a typical circular 
velocity of ar ~ |Ai;|/2, then ar — 1100 km s^^, which does 
not exceed the sound velocity of the Coma cluster, which is 
1500 km s~^ for a temperature of 8 keV. The average redshift 
in the center region is as follows: 0.0227± 0.0011 in the north- 
west region (ceU numbers 8, 9, 12, and 13), 0.0214 ± 0.0013 
in the south-west region (cell numbers 10, 11, 14, and 15), 
0.0236 ± 0.0012 in the north-east region (cell numbers 0, 1, 4, 
and 5), and 0.0222 ±0.0012 in the south-east region (cell num- 
bers 2, 3, 6, and 7). We also show the averaged redshift over 
FOV in the upper panel in table 4. These results are consis- 
tent with each other in the center region, and with the observed 
redshifts in the 14' offset region within statistical and system- 
atic errors. The difference in the ICM velocity between the 
center and 30' offset regions is derived as Av = 90± 1350 ± 818 
km s~^, and the difference between the center and 34' offset 
regions is Av — —1050 ± 3300 ± 818 km s^^. These are also 
consistent with each other. Components of the velocity in the 
plane of the sky are not measured. If we assume that three 
components are similar, then the limit on the radial velocity 
translates into a limit of \/3 x 290 ~ 500 km s^^, which does 
not exceed the speed of sound and is consistent with systematic 
errors. Therefore, there were no significant velocity variations 
on the 10' scale within central regions (center and 14', 30', 
and 34' offset regions) from our analysis. In the NGC 4839 
subcluster region (44' offset region), the observed redshift is 
0.0269±0.0033, which corresponds to 1140±990 km s"^ as the 
ICM bulk velocity without systematic errors. There also seems 
to be no significant difference between the core of the Coma 
cluster and the NGC 4839 subcluster within the 90% confi- 
dence level. 

As mentioned in subsection 4.5.1, some sensors have large 



uncertainties in the determination of the line centroid energy. 
On XIS3 in the center region, the Mn Ka line centroid energy 
of segment A was different from others. Similarly, the segment 
D of XISl in the NGC 4839 region was also different. We av- 
eraged the spectrum over FOV in (i) all detectors, and (ii) all 
detectors except those with different Mn Ka line centroid en- 
ergies to check the spectral shape and measured redshift in the 
center and NGC 4839 regions (44' offset region). The results 
are shown in the lower section in table 4. For the center region, 
the derived redshifts were fairly consistent including XIS3 or 
not. On the other hand, for the NGC 4839 subcluster region, 
the resultant redshift changed slightly from 0.0268 to 0.0287, 
although this was also consistent within statistical and sys- 
tematic errors. As shown in figure 11, the peak of the He-like 
Fe-Ka line of XISl is more shifted than that of the other two 
sensors. The measured redshifts in the 44' offset region were 

0. 0271±0.0023, 0.0303±0.0051, and 0.0255±0.0027 for XISO, 

1, and 3, respectively. Because the measured redshift of the 
averaged spectra between XISO and 3 was consistent with the 
redshift from simultaneous spectral fits, the value 0.0269 would 
be plausible. 

In conclusion, there was no significant difference in the 
ICM velocity between the core of the Coma cluster and the 
NGC 4839 subcluster within 90% confidence. 

4-6. Fe abundance of ICM 

To measure the abundance using the strong Fe-K lines, we 
employed the results of spectral fitting in the 5.0-8.0 keV en- 
ergy range for the Fe abundance. Because of the good statis- 
tics, the spectra in the lower energy band primarily determine 
the fitted temperatures. Although the spectra are well fitted 
with a single-temperature model, a small discrepancy between 
the model and data around the Fe-K line can yield a large 
systematic uncertainty in the derived Fe abundances. Within 
the 34' region, temperatures and abundances derived from the 
spectral fitting in the 5.0-8.0 keV energy range agree well with 
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Table 4. Results of the redshift measurement with the single temperature (APEC) model in the energy range of 5.0-8.0 keV. The 
upper section shows the averaged redshift for each region with all XIS sensors simultaneously. The lower section shows the redshift in 
each field with the averaged spectral fits. XIS detectors that are included are shown in parentheses. Systematic errors of ±818 km s~^ 
in the X-ray (cz) were not included in this table. Optical rcdshifts are referred in CoUess & Dunn (1996). 



field 


X-ray (z) 


Optical (z) 


X-ray (cz) 
km s 


Optical (cz) 


all pointings 


0.0227±0.0005 


0.02307 


6810±150 


6917 


center 


0.0227±0.0005 


0.02286 


6810±150 


6853 


14' offset 


0.0230±0.001 


0.02286 


6900±300 


6853 


center and 14' offset 


0.0227±0.0005 


0.02286 


6810±150 


6853 


NGC 4839 


0.0269±0.0033 


0.02448 


8070±990 


7339 


center (XIS0123) 
center (XIS012) 
NGC 4839 (XIS013) 
NGC 4839 (XIS03) 


0.0226±0.0005 
0.0226±0.0005 
0.0287±0.0033 
0.0268±0.0033 


0.02286 
0.02286 
0.02448 
0.02448 


6775±150 
6775±150 
8604±990 
8034±990 


6853 
6853 
7339 
7339 
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Fig. 13. Left: Radial Fe abundance profile with the single temperature (APEC) model in the 5.0—8.0 keV energy range measured 
from the position of the X-ray peak of the Coma cluster. Color and symbol notations are the same as in figure 6. The gray colors 
show the results from XMM-Newton (Matsushita 2011). Right: Observed Fe abundance map of the Coma cluster in solar units. 
The magenta circle corresponds to the X-ray peak. The Fe abundance map is overlaid on X-ray contours on a linear scale in the 
0.4-2.4 keV energy range with the ROSAT-AU-Sky-Survey. 



those from 1.0-8.0 keV energy ranges within 5%. In the 44' 
and 60' offset regions, the Fe abundances from the 5.0-8.0 keV 
energy range are smaller than those from the 1.0-8.0 keV range 
by ~0.1 solar, values for the limited energy range of 5.0- 
8.0 keV from the former fitting are 10% smaller than those 
from the latter, although the single-temperature model repro- 
duces the observed spectra at 1.0-8.0 keV well. Therefore, we 
adopted the Fe abundance derived from the spectral fitting in 
the 5.0-8.0 keV energy range. Figure 13 shows a radial profile 
and map of the Fe abundance of ICM with the single temper- 
ature (APEC) model in the 5.0-8.0 keV energy range. Within 
30' from the X-ray peak, Fe abundances are almost constant 
at ~0.4 solar. The derived Fe abundances are 0.41±0.01 solar 
in the center region, and 0.38±0.02 solar in the 14' offset re- 
gion. The Fe abundances in the 30' and 34' offset regions are 
slightly lower, ~0.3 solar, although it has larger uncertainties. 
These values agree with recent XMM results (Arnaud et al. 
2001; Matsushita 2011) within 1.2 Mpc. 

The Fe abundance in the NGC 4839 subcluster region is 



comparable to that in the central region, and higher than that 
in the 30' and 34' offset regions. The average value in the sub- 
cluster region is 0.33±0.09 solar. This value is also consistent 
with that in Neumann et al. (2001), although the error bars in 
their results are much larger than those in our results. 

5. Discussion 

5.1. Temperature measurements with Suzaku 

X-ray observations of clusters of galaxies can play a ma- 
jor role in determining cosmological parameters, because of 
the importance of measuring gravitational mass. The ratio of 
He-like Kq to H-like Ka lines of Fe is a very steep function of 
temperature, and has a different temperature dependence from 
that of the continuum. The continuum emission is affected by 
the systematic uncertainties in the response and non-thermal 
emission, while the line ratio does not. Thus, by comparing 
the fitted temperatures between these two methods, the effects 
of the systematic uncertainty in the response and/or of the 
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non-thermal emission is increased. Using the APEC plasma 
code, the temperatures derived from the spectral fitting of the 
1.0-8.0 keV energy range and the line ratio of Fe agree within 
a few percent when we averaged all regions. The systematic 
uncertainties in temperatures from the line ratio are in the 
level of the systematic difference between those from APEC 
and MEKAL, which is about 3%. The consistency of observed 
temperatures from the continuum and those from line ratios 
supports the accuracy of Suzaku measurements of plasma tem- 
perature. In contrast, systematic differences of ~ 10% in the 
cluster temperature among Chandra, PN, and MOS were re- 
ported in relatively high-temperature clusters (Snowden et al. 
2008). Even using the same instrument, the derived temper- 
atures changed systematically by ~ 10% between calibrations 
(Reese et al. 2010; Matsushita 2011). 

A systematic uncertainty of several percent in response ma- 
trices of a detector can also cause systematic uncertainties in 
the temperature structure. Matsushita (2011) found that the 
temperature structures of ICM in the Coma cluster derived 
from multi-temperature fits with PN and MOS spectra were 
different, although the two detectors gave consistent tempera- 
tures within several percent. When they allowed higher tem- 
perature components over 10 keV in hot clusters including the 
Coma cluster, using MOS data, Fe abundances from the He- 
like line increased by several tens of percent, owing to artificial 
detection of a hot component, while those derived by PN were 
unchanged. However, with Suzaku, we found that the single- 
temperature model fits the spectra of the Coma cluster very 
well, and we do not need temperature components above 10 
keV or below 2 keV to fit the spectra. The multi-temperature 
components observed by XMM were not needed. 

In the outer regions of clusters at low surface brightness, 
uncertainties in the background of Chandra and XMM also 
lead to systematic uncertainties in ICM temperatures. Some 
analyses found that ICM temperatures decreased outward by 
half at 0.6ri8o (e.g., Vikhlinin et al. 2005; Pratt et al. 2007), 
while others found flatter profiles (e.g., Arnaud et al. 2005). 
Systematic uncertainties in the derived temperatures due to 
uncertainties in the background with Suzaku observations are 
much smaller, because Suzaku XIS has a much lower and more 
stable background. Also, previous XMM and Chandra results 
of temperatures in the outer regions have statistical errors more 
than 20%, while those of Suzaku have errors less than 10%. 

5.2. Effect on the search of a hard component 

We need precise measurements of ICM temperatures to 
search for a non-thermal hard component as an excess of the 
thermal component. On the basis of the ICM temperature 
map derived with XMM-PN, Wik et al. (2009) found that the 
Suzaku HXD-PIN spectrum was described with the thermal 
emission from ICM, and there was no significant evidence of 
excess hard emission. The temperatures of the Coma cluster 
derived from Suzaku-XIS are consistent within error bars with 
those from the PN detector: the weighted average of temper- 
atures with statistical errors within the center field is close to 
that derived with XMM-PN within a few percent. The temper- 
atures of offset fields are consistent with those from XMM-PN 
within error bars. Thus, our results support the non-detection 
of non-thermal hard emissions by Wik et al. (2009) 

Ota et al. (2008) reported the existence of hot gas at 25.3l|'5 
keV in RX J1347.5-1145. They concluded that the gas prop- 
erties can be explained by a fairly recent (within the last 0.5 
Gyr) collision of two massive (5x10^* M© ) clusters with bullet- 
like high velocity (An ~ 4500 km s~^). In contrast, Wik et al. 



(2009) found that there was no extremely high temperature 
gas in the Coma cluster in the HXD-PIN spectrum up to 70 
keV. On the basis of the consistency of derived temperatures 
from the continuum and the line ratio and the results from the 
spectral analysis of HXD-PIN, we conclude that the amount 
of very high temperature gas is rather small. From the spatial 
variation of the derived temperatures, there is no region with 
extremely hot temperature: the derived temperatures range 
from 7.0 to 9.0 keV within the center, 14' offset, and 34' offset 
regions. 

The Coma cluster has a cluster-wide synchrotron radio halo. 
If the electrons in the cluster were accelerated by a very strong 
shock that occurred owing to a recent merger, we might expect 
the cluster to contain more very hot gas, although it is difficult 
to make this argument quantitative. In any case, in the Coma 
cluster, relativistic electrons producing the radio halo were ap- 
parently generated without producing a large amount of very 
hot gas. An example of a process which might do this is the 
turbulent re-acceleration model (e.g., Brunetti et al. 2001). 

5.3. Temperature structure of the Coma cluster 

Within the center and 14' offset regions the spatial variation 
of derived temperatures is relatively small, from 7 keV to 9 
keV. There are no extremely hot and cool regions. We did not 
detect any temperature jump corresponding to a shock in the 
offset regions. Within each region for the spectral analysis, the 
single-temperature model fit the spectrum. Therefore, in the 
center of the Coma cluster, ICM is in a nearly relaxed state. 

Although the variation was small, there was some asym- 
metric temperature structure, like a hot region northwest of 
the center. If the cluster had been an isolated system col- 
lapsing from self-gravity, the temperature structures would be 
symmetric. The observed asymmetry in the temperature dis- 
tribution indicates that external effects or interactions, such 
as merging, have occurred. Numerical simulations of the evo- 
lution of the cluster (e.g., Roettiger et al. 1996; Norman & 
Bryan 1999) show that nonaxisymmetric structures in the tem- 
perature distribution are erased several Gyr after each merger. 
The presence of asymmetric temperature structures suggest 
that the central region of the Coma cluster has experienced a 
merger, and the lack of hard component shows that its core 
have relaxed in the last several Gyr. 

The observed cool region in the southeast of the center re- 
gion coincides with the filamentary structure originating near 
NGC 4911 and NGC 4921 (Vikhlinin et al. 1997; Donnelly et 
al. 1999). Arnaud et al. (2001) reported that the temperature 
of this cold area ranges from 5 to 7 keV. This region mostly 
contain the stripped subcluster gas, which is thought to be rela- 
tively cool in the initial state. The numerical simulations (e.g., 
Schindler & Miiller 1993) suggest that the gas would maintain 
its temperature for several Gyr. Ishizaka & Mineshige (1996) 
showed that even a single subcluster collision with supersonic 
velocity creates both hot and cool regions in the main cluster. 
Although heating due to collisions occurs, cooling would also 
occur through the adiabatic expansion of the stripped gas. 

5.4. Searching for gas bulk motions 

The measurements of bulk motions in ICM are also impor- 
tant to derive the gravitational mass of clusters of galaxies. 
We constrained the upper limit of ICM bulk motions in the 
center and 14' offset regions on the scale of 117 kpc (4.'5) and 
234 kpc (9'), respectively, which correspond to the scales of 
regions for spectral accumulation (see section 2 and subsection 
4.5.3). The upper limit of gas bulk motions is \Av\ < 1100 
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km s~^. As summarized in table 4, considering the system- 
atic error of 18 eV, or 818 km s^^ with a 90% confidence level 
(Section 4.5.1), there are no significant bulk motions in the 
center and 14' offset fields. Compared with the sound velocity 
of the Coma cluster, s = {'ykT/ ump)^^^ ~ 1500 km s~^ for a 
temperature of kT — 8 keV, the bulk-motion in the central 
region of the Coma cluster does not exceed the sound velocity. 

Our results do not conflict with the discovery of pressure 
fluctuations in the center of the Coma cluster by Schuecker et 
al. (2004) which indicates that at least 10% of the ICM pressure 
is in turbulent form. We have not constrained the turbulence, 
and the scale of the fluctuation ranges between 29 kpc and 64 
kpc, which is smaller than our scale of spectral analysis. 

The velocity difference between the main cluster and the 
NGC 4839 subcluster from optical observations is 486 km s~^. 
Although the best-fit line-of-sight velocity of the subcluster is 
higher than that of the main cluster by ~ 1000 km s~^, the 
difference is comparable to statistical and systematic errors. 
Because the temperature of the region is about 5 keV, the 
sound velocity is ~ 1150 km s~^. Considering the error bars, 
we conclude that the relative line-of-sight gas velocity in the 
NGC 4839 subgroup region does not greatly exceed the sound 
velocity. 

Numerical simulations indicate that a major merger between 
clusters raises the temperature of ICM and induces a bulk ve- 
locity of the order of 1000 km s^^ (e.g., Roettiger et al. 1996; 
Norman & Bryan 1999). In the late phase of mergers, turbu- 
lent motions develop in a simulation by Takizawa (2005). The 
non-detection of bulk-motions in ICM of the Coma cluster also 
indicates that the central region of the cluster is in a somewhat 
relaxed state. If so, the non-thermal electrons relevant to the 
radio halo of the Coma cluster may have been accelerated by 
the intracluster turbulence rather than shocks. In conclusion, 
the line-of-sight gas velocity of the Coma cluster in the ob- 
served region does not greatly exceed the sound velocity, and 
the central region of the Coma cluster is in a nearly relaxed 
state. This means that the assumption of hydrostatic equilib- 
rium is approximately valid in calculating the cluster mass of 
the Coma cluster. 

5.5. Fe Abundance map 

The Fe abundance is nearly constant at ~ 0.4 solar within 
a clustocentric radius of 0.5 Mpc, and beyond that distance 
it decreases with the radius to 0.2-0.3 solar. This profile is 
consistent with those measured by XMM (Arnaud et al. 2001; 
Matsushita 2011). Suzaku has a lower level of background, 
and a slightly better energy resolution: therefore, systematic 
uncertainties in Fe abundance are considered to be smaller. 

Differences in abundance profiles of ICM in central regions 
of clusters with and without cool cores have been reported 
(Fukazawa et al. 2000, De Grandi & Molendi 2001, Matsushita 
2011). Matsushita (2011) showed an abundance difference be- 
tween the two type of clusters within O.lrigo, while at 0.1- 
0.3ri8o, the average values of Fe abundances are both ~ 0.4-0.5 
solar, and are consistent with each other. 

The Coma cluster is thought to have experienced a major 
merger in the recent past and subsequently, the central region 
came to a somewhat relaxed state. During a cluster merger, 
the mixing of ICM could destroy any central Fe peak. In the 
first stage of merging, the Fe peak and cool core remain in- 
tact as observed in some ongoing merging clusters (Molendi et 
al. 2000; Sun et al. 2002). For example, Abell 2256, an ongo- 
ing merging cluster, has an abundance gradient in the direc- 
tion opposite of the merging direction owing to the remixing 



of the gas. Numerical simulation also support this scenario. 
Schindler et al. (2005) performed hydrodynamical simulations 
with two different merger types of clusters: Cluster 1 has only 
small merger events, while Cluster 2 undergoes a major merger 
(mass ratio 1:3). In Cluster 1, the abundance is more than ~ 
0.4 solar within a radius of 0.1 Mpc, and decreases from ~ 0.3 
solar to ~ 0.1 solar with a decrease in radius. The situation 
is different in Cluster 2. In the Coma cluster, the Fe abun- 
dance is ~ 0.4 solar and homogeneous within a radius of 0.5 
Mpc. This value of the Fe abundance is similar to those at 
0.2-0.3ri8o in other clusters as shown in figure 13. Therefore, 
tn the Coma cluster, the merger have destroyed the cool core, 
and mixed the ICM completely at least within 0.5 Mpc. 

6. Summary &: Conclusion 

The Coma cluster which was observed with the X-ray 
Imaging Spectrometer (XIS) onboard Suzaku, was analyzed 
by the X-ray satellite in six pointings, centered on the X-ray 
peak and offset by 14' west, 30', 34', 44', and 60'. Because of 
its low background, Suzaku is the most sensitive X-ray satellite 
for observing Ka lines of Fe in the intracluster medium. After 
obtaining accurate measurements of Fe lines, we studied the 
temperature structure of the intracluster medium, searched for 
possible bulk-motions, and measured the Fe abundance in the 
cluster. 

The spectra of each extracted region were well fitted by the 
single-temperature APEC model, and the two- or three- tem- 
perature APEC model did not improve x^- The temperatures 
derived from the observed ratios of Ka lines of H-like and He- 
like Fe agree well with those from spectral fittings with the 
single-temperature APEC model. Because this line ratio is 
a strong function of plasma temperature, the observed con- 
sistency supports the accuracy of temperature measurements 
with Suzaku, and constrains the temperature structure of ICM. 

The observed values of the central energy of the He-like Fe 
line of the center, 14', 30' and 34' offset regions are constant 
within 500 km s^^, which corresponds to the calibration error. 
Because relative bulk velocities in the Coma cluster are smaller 
than the sound velocity of the intracluster medium, 1500 km 
s we can verify the derived total mass in a cluster on the 
basis of the hydrostatic ICM equilibrium inside the 34' offset 
region. Significant bulk velocities were also not found in the 44' 
offset region, which corresponds to the NGC 4839 subcluster. 

The results on the temperature and velocity structure sug- 
gest that the core of the Coma cluster is in a fairly relaxed 
state. This is consistent with models in which the non-thermal 
electrons relevant to the radio halo are accelerated by the in- 
tracluster turbulence rather than large-scale shocks. This is 
also consistent with the fact that numerical simulations show 
that the turbulence motion is developed in the late phase of 
mergers. 

The observed Fe abundance of the intracluster medium is 
almost constant at 0.4 solar inside the 34' offset region, and 
decreases with radius. The central abundance is slightly lower 
than that of other clusters or groups. This indicates that 
central regions of the gas were mixed well during the past 
merger growth of the cluster. 
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